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Abstract 

This paper presents two new constructions related to singular solutions of polynomial 
systems. The first is a new deflation method for an isolated singular root. This construction 
uses a single linear differential form defined from the Jacobian matrix of the input, and defines 
the deflated system by applying this differential form to the original system. The advantages 
of this new deflation is that it does not introduce new variables and the increase in the 
number of equations is linear instead of the quadratic increase of previous methods. The 
second construction gives the coefficients of the so-called inverse system or dual basis, which 
defines the multiplicity structure at the singular root. We present a system of equations in the 
original variables plus a relatively small number of new variables. We show that the roots of 
this new system include the original singular root but now with multiplicity one, and the new 
variables uniquely determine the multiplicity structure. Both constructions are “exact” in that 
they permit one to treat all conjugate roots simultaneously and can be used in certification 
procedures for singular roots and their multiplicity structure with respect to an exact rational 
polynomial system. 


1 Introduction 

Our motivation for this work is twofold. On one hand, in a recent paper If], two of the co-authors 
of the present paper studied a certification method for approximate roots of exact overdetermined 
and singular polynomial systems, and wanted to extend the method to certify the multiplicity 
structure at the root as well. Since all these problems are ill-posed, in pQ a hybrid symbolic- 
numeric approach was proposed, that included the exact computation of a square polynomial 
system that had the original root with multiplicity one. In certifying singular roots, this exact 
square system was obtain from a deflation technique that added subdeterminants of the Jacobian 
matrix to the system iteratively. However, the multiplicity structure is destroyed by this deflation 
technique, that is why it remained an open question how to certify the multiplicity structure of 
singular roots of exact polynomial systems. 

Our second motivation was to find a method that simultaneously refines the accuracy of a singular 
root and the parameters describing the multiplicity structure at the root. In all previous numerical 
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approaches that approximate these parameters, they apply numerical linear algebra to solve a linear 
system with coefficients depending on the approximation of the coordinates of the singular root. 
Thus the local convergence rate of the parameters was slowed from the quadratic convergence of 
Newton’s iteration applied to the singular roots. We were interested if the parameters describing 
the multiplicity structure can be simultaniously approximated with the coordinates of the singular 
root using Newton’s iteration. 

In the present paper we first give a new improved version of the deflation method that can be 
used in the certification algorithm of [T], reducing the number of added equations at each deflation 
iteration from quadratic to linear. We prove that applying a single linear differential form to the 
input system, corresponding to a generic kernel element of the Jacobian matrix, already reduced 
both the multiplicity and the depth of the singular root. Secondly, we give a description of the 
multiplicity structure using a polynomial number of parameters, and express these parameters 
together with the coordinates of the singular point as the roots of a multivariate polynomial system. 
We prove that this new polynomial system has a root corresponding to the singular root but now 
with multiplicity one, and the new added coordinates describe the multiplicity structure. Thus this 
second approach completely deflates the system in one step. The number of equations and variables 
in the second construction depends polynomially on the number of variables and equations of the 
input system and the multiplicity of the singular root. Both constructions are exact in the sense 
that approximations of the coordinates of the singular point are only used to detect numerically 
non-singular submatrices, and not in the coefficients of the constructed polynomial systems. 
Related work. 

The treatment of singular roots is a critical issue for numerical analysis and there is a huge literature 
on methods which transform the problem into a new one for which Newton-type methods converge 
quadratically to the root. 

Deflation techniques which add new equations in order to reduce the multiplicity have already 
been considered in 1261 . [25] : By triangulating the Jacobian matrix at the (approximate) root, new 
minors of the polynomial Jacobian matrix are added to the initial system in order to reduce the 
multiplicity of the singular solution. 

A similar approach is used in [lOj and [2, where a maximal invertible block of the Jacobian matrix 
at the (approximate) root is computed and minors of the polynomial Jacobian matrix are added 
to the initial system. In [8] , an additional step is considered where the first derivatives of the input 
polynomials are added when the Jacobian matrix at the root vanishes. 

These constructions are repeated until a system with a simple root is obtained. 

In these methods, at each step, the number of added equations is (n — r) x (m — r), where n is 
number of variables, m is the number of equations and r is the rank of the Jacobian at the root. 
In |12] , a triangular presentation of the ideal in a good position and derivations with respect to the 
leading variables are used to iteratively reduce the multiplicity. This process is applied for p-adic 
lifting with exact computation. 

In other approaches, new variables and new equations are introduced simultaneously. In m, 
new variables are introduced to describe some perturbations of the initial equations and some 
differentials which vanish at the singular points. This approach is also used in [ITS], where it is 
shown that this iterated deflation process yields a system with a simple root. 

In [20] , perturbation variables are also introduced in relation with the inverse system of the singular 
point to obtain directly a deflated system with a simple root. The perturbation is constructed from 
a monomial basis of the local algebra at the multiple root. 

In [T3] Hi] , only variables for the differentials of the initial system are introduced. The analysis 
of this deflation is improved in [5] , where it is shown that the number of steps is bounded by the 
order of the inverse system. This type of deflation is also used in EH, for the special case where 
the Jacobian matrix at the multiple root has rank n — 1 (case of breath one). 

In these methods, at each step, the number of variables is at least doubled and new equations are 
introduced, which are linear in these new variables. 

The mentioned deflation techniques usually breaks the structure of the local ring at the singular 
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point. The first method to compute the inverse system describing this structure is due to F.S. 
Macaulay pj3] and known as the dialytic method. More recent algorithms for the construction of 
inverse systems are described e.g. in |2lj . reducing the size of the intermediate linear systems (and 
exploited in [28]). In [32], the dialytic method is used, and they analyze the relationship of deflation 
some methods to the inverse system. It had been further improved in 1231 and more recently in 
[20] . using an integration method. This technique reduces significantly the cost of computing 
the inverse system, since it relies on the solution of linear system related to the inverse system 
truncated in some degree and not on the number of monomials in this degree. Multiplication 
matrices corresponding to systems with singular roots were studied in |221 4J. 

The computation of inverse systems has been used to approximate a multiple root. In 122], a 
minimization approach is used to reduce the value of the equations and their derivatives at the 
approximate root, assuming a basis of the inverse system is known. In [30] . the inverse system 
is constructed via Macaulay’s method; tables of multiplications are deduced and their eigenvalues 
are used to improve the approximated root. They show that the convergence is quadratic at the 
multiple root. In [16, they show that in the breadth one case the parameters needed to describe 
the inverse system is small, and use it to compute the singular roots in [15] . In [20], the inverse 
system is used to transform the singular root into a simple root of an augmented system. 
Contributions. 

We propose a new deflation method for polynomial systems with isolated singular points, which 
does introduce new parameters. At each step, a single differential of the system is considered based 
on the analysis of the Jacobian at the singular point. A linear number of new equations is added 
instead of the quadratic increases of the previous deflations. The deflated system does not involved 
any approximate coefficients and can therefore be used in certification methods as in 1]. 

To approximate efficiently both the singular point and its inverse system, we propose a new defla¬ 
tion, which involves a small number of new variables compared to other approaches which rely on 
Macaulay matrices. It is based on a new characterization of the isolated singular point together 
with its multiplicity structure. The deflated polynomial system exploits the nilpotent and commu¬ 
tation properties of the multiplication matrices in the local algebra of the singular point. We prove 
that it has a simple root which yields the root and the coefficients of the inverse system at this 
singular point. Due to the upper triangular form of the multiplication matrices in a convenient 
basis of the local algebra, the number of new parameters introduced in this deflation is less than 
i n(S — 1)J where n is the number of variables and J the multiplicity of the singular point. The 
parameters involved in the deflated system are determined from the analysis of an approximation of 
the singular point. Nevertheless, the deflated system does not involve any approximate coefficients 
and thus it can also be used in certification techniques as |T]. 

In this paper we present two new constructions. The first one is a new deflation method for a 
system of polynomials with an isolated singular root. The new construction uses a single linear 
differential form defined from the Jacobian matrix of the input, and defines the deflated system 
by applying this differential form to the original system. We prove that the resulting deflated 
system has strictly lower multiplicity and depth at the singular point than the original one. The 
advantage of this new deflation is that it does not introduce new variables, and the increase in 
the number of equations is linear, instead of the quadratic increase of previous deflation methods. 
The second construction gives the coefficients of the so called inverse system or dual basis, which 
defines the multiplicity structure at the singular root. The novelty of our construction is that 
we show that the nilpotent and commutation properties of the multiplication matrices define 
smoothly the singular points and its inverse system. We give a system of equations in the original 
variables plus a relatively small number of new variables, and prove that the roots of this new 
system correspond to the original multiple roots but now with multiplicity one, and they uniquely 
determine the multiplicity structure. The number of unknowns used to describe the multiplicity 
structure is significantly smaller, compared to the direct computation of the dual bases from the 
so called Macaulay matrices. Both constructions are “exact” in the sense that approximations of 
the coordinates of the singular point are only used to detect numerically non-singular submatrices, 
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and not in the rest of the construction. Thus these constructions would allow to treat all conjugate 
roots simultaneously, as well as to apply these constructions in the certification of the singular 
roots and the multiplicity structure of an exact rational polynomial system. 


2 Preliminaries 

Let f := (/i,..., /at) G K[x] w with x = (xi ,..., x n ) for some IcC held. Let £ = (£ 1 ,..., £„) G 
C™ be an isolated multiple root of f. Let / = (/i,..., /at), be the maximal ideal at £ and Q be 
the primary component of / at £ so that yfQ = . 


Consider the ring of power series C[[c^]] := C[[c?i,£,..., c^]] and we use the notation for (3 = 


We identify C [[#$]] with the dual space C[x]* by considering d^ as derivations and evaluations at 
£, dehned by 


af(p) := d P {p)\i := 




dXi 1 ■ ■ ■ dXn 


-(C) 


for p G C[x]. 


(1) 


Hereafter, the derivations “at x” will be denoted d 13 instead of d The derivation with respect to 
the variable d t is denoted dg i {i = 1,..., n). Note that 


1 


0|((x-0“) 


1 it a = (3 
0 otherwise ’ 


where we use the notation ~ j. 1 .^ , . 

For p G C[x] and A G C[[c^]] = C[x]*, let 

p ■ A : q K > A (pq). 


We check that p = (xi — &) acts as a derivation on C[[0^]]: 

(Xi -&) ■ d | = dg i( {d^) 

For an ideal I C C[x], let I 1 = {A G C[[c?{]] | Vp G /, A(p) = 0}. The vector space I 1 - is naturally 
identified with the dual space of C[x]/J. We check that I 3 - is a vector subspace of C[[0^]], which 
is stable by the derivations dg i ( . 

Lemma 2.1. If Q is a m c-primary component of I, then Q 3 - = I 3 - (1 C[0|]. 

This lemma shows that to compute Q 1 -, it suffices to compute all polynomials of C[9j] which are 
in I- 1 . Let us denote this set S = I 3 - C I C[9$]. It is a vector space stable under the derivations 
dg i( . Its dimension is the dimension of Q 3 - or C[x]/<2, that is the multiplicity of £, denote it by 
<5j(J), or simply by S if £ and I is clear from the context. 

For an element A(d^) G C[9$] we define the order ord(A) to be the maximal \(3\ such that d | 
appears in A(dc) with non-zero coefficient. 

For t G N, let St be the elements of S of order < t. As S is of dimension d. there exists a smallest 
t > 0 such that St+i = St- Let us call this smallest t, the nil-index of St and denote it by o^(I), 
or simply by o. As S is stable by the derivations dg i ( , we easily check that for t > 0^(1), St = S 
and that o^(I) is the maximal degree of the elements in S. 
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3 Deflation using first differentials 


To improve the numerical approximation of a root, one usually applies a Newton-type method to 
converge quadratically from a nearby solution to the root of the system, provided it is simple. In 
the case of multiple roots, deflation techniques are employed to transform the system into another 
one which has an equivalent root with a smaller multiplicity or even with multiplicity one. 

We describe here a construction, using differentials of order one, which leads to a system with a 
simple root. This construction improves the constructions in mm since no new variables are 
added. It also improves the constructions presented in pJD| and the “kerneling” method of El by 
adding a smaller number of equations at each deflation step. Note that, in [Sj, there are smart 
preprocessing and postprocessing steps which could be utilized in combination with our method. 
In the preprocessor, one adds directly partial derivatives of polynomials which are zero at the root. 
The postprocessor extracts a square subsystem of the completely deflated system for which the 
Jacobian has full rank at the root. 

Consider the Jacobian matrix Jf (x) = [d :/ fi(x)\ of the initial system f. By reordering properly the 
rows and columns (i.e., polynomials and variables), it can be put in the form 


J f (x) 


A(x) B(x) 
C(x) D(x) 


( 2 ) 


where A(x) is an r x r matrix with r = rankJf(£) = rankA(£). 
Suppose that B(x) is anrxc matrix. The c columns 


det(A(x)) 


-A- 1 (x)B(x) 

Id 


(for r = 0 this is the identity matrix) yield the c elements 

n n 

A? = £ A, j (x)<9j , ..., A* = £ X c j (x)<9j . 

*=l i=i 

Their coefficients A^-(x) € K[x] are polynomial in the variables x. Evaluated at x = £, they 
generate the kernel of Jf(£) and form a basis of 

Definition 3.1. The family Df = {Af ,..., A*} is the formal inverse system of order 1 at £. For 
* = {*!,..., ik} C {1,..., c} with |*| 7 ^ 0, the i-deflated system of order 1 of f is 

{f,A*(f),...,A*(f)}. 


By construction, for i = 1,..., c, 

n 

A ?( f ) = 5Z^( f ) Ai d( x ) = det(A(x))J f (x)[A M (x)] 
j=i 

has n — c zero entries. Thus, the number of non-trivial new equations added in the *-deflated 
system is |*| • (N — n + c). The construction depends on the choice of the invertible block A(£) in 
Jf (£). By a linear invertible transformation of the initial system and by computing a *-deflated 
system, one obtains a deflated system constructed from any |*| linearly independent elements of 
the kernel of Jf(£)- 

Example 3.2. Consider the multiplicity 2 root £ = (0,0) for the system /i(x) = x\ + x\ and 
/ 2 (x) = xf + x\. Then, 

1 2*2 ' 

2*i 2*2 


Jf(x) = 


A(x) 

B(x) 


C(x) 

£>(x) 
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The corresponding vector [—2x2 1] T yields the element 

A* = -2x 2 d! + d 2 . 

Since A*(/i) = 0, the {l}-deflated system of order 1 of f is 

{x\+x\, xf + x%, —4 xiX 2 +2 x 2 } 
which has a multiplicity 1 root at £. 

We use the following to analyze this deflation procedure. 

Lemma 3.3 (Leibniz rule). For a,b £ K[»], 

d a (ab) = J^{a)d^d a ){b). 

p ' 

Proposition 3.4. Let r be the rank of Jf(£). Assume that r < n. Let i C {l,...,n} with 
0 < \i\ < n — r and be the *-deflated system of order 1 of f. Then, yfU)) > 1 and 
oe(fU)) < oj(f), which also implies that <5c(f (1 l) < de(f). 

Proof. By construction, for i £ i, the polynomials Af(f) vanish at f, so that > 1. By 

hypothesis, the Jacobian of f is not injective yielding oj(f) > 0. Let be the inverse system 
of fA) at f. Since (f^) D (f), we have c In particular, for any non-zero element 

A £ K[c>^] and i £ i, A(f) = 0 and A(Af (f)) = 0. 

Using Leibniz rule, for any p £ K[x], we have 




A ( X ) a 3 (p) 


( A <.i( x )) d f £ ( A ) a i,«(p) 

/3€N n j = l 


53 Yf 4 a f( A ci( x )) a 3.« d ?,( A )(p) 

pern j=i 

^ A ii@ df (A)(p) 
pe 


where 

n ^ ^ 

A= y K,j,pdj£ £ K[0 e ] and A* = —df (Xi,j(x)) £ K. 

j = 1 

The term A,.o is 33j=i i which has degree 1 in cL since [A»,j(£)] is a non-zero element of 

ker J f (£). For simplicity, let fa (A) := A i)/3 df(A). 

For any A £ C[0^], we have 

dg J( (fa(A)) = y Ki,pdo(A) + Ai^d§(d dji (A)) 

peN n 

= y ^i,jpdg{A) + fa(d dji (A)). 


Moreover, if A £ then by definition 0j(A)(f) = 0. Since @ and S>A) a re both stable by 

derivation, it follows that VA £ & 1 ), dg ji (fa( A)) £ fiA 1 ) + fa(0&^). Since ( 2 ) { A') C Q>, we know 
+ fa{2>^) is stable by derivation. For any element A of + fa{S>^), A(f) = 0. We deduce 
that Q> + fa(S>^) = Q>. Consequently, the order of the elements in fa(@^) is at most oc(f). The 
statement follows since fa increases the order by 1, therefore ot(f^^) < oj(f). □ 
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We consider now a sequence of deflations of the system f. Let f 1 W be the ii-deflated system of f. We 
construct inductively f( fc+1 ) as the ifc+i-deflated system of for some choices of ij C {1, ..., n}. 

Proposition 3.5. There exists k < o^(f) such that £ is a simple root of fl fe ). 

Proof. By Proposition 13.41 (5j(f( fc )) > \ and is strictly decreasing with k until it reaches 

the value 0. Therefore, there exists k < oj(/) such that o^(f( fe l) = 0 and > 1. This implies 

that £ is a simple root of fw. □ 

To minimize the number of equations added at each deflation step, we take |i| = 1. Then, the 
number of non-trivial new equations added at each step is at most N — n + c. 

We described this approach using first order differentials arising from the Jacobian, but this can 
be easily extended to use higher order differentials. 

4 The multiplicity structure 

Before describing our results, we start this section by recalling the definition of orthogonal primal- 
dual pairs of bases for the space C[x]/Q and its dual. The following is a definition/lemma: 

Lemma 4.1 (Orthogonal primal-dual basis pair). Let f, f, Q, 3>, 8 = J{(f) and o = o^(f) he 
as above. Then there exists a primal-dual basis pair of the local ring C[x]/Q with the following 
properties: 

• The primal basis of the local ring C[x]/Q has the form 

B := {(x - 0 Q1 , (x - £)“ 2 , • • •, (x - £)“*}• (3) 


We can assume that ol\ = 0 and that the monomials in B are connected to 1 (c.f. \24\ I). 
Define the set of exponents in B 


E := {ai,... ,as}. 


• There is a unique dual basis Ac^ orthogonal to B, i.e. the elements of A are given in the 
following form: 


Ao 

Ai 


=u 

-Vr + 

a\\ * 


5Z d t 

m <o 


Aa_i 


—df 

cts- 4 


^ 1 u as,/3 9^ i 
\P\<o 

(IgE 


(4) 


• We have 0 = ord(Ao) < • • • < ord(Aj_i), and for all 0 < t < o we have 

@t = span{Aj : ord(Aj) < t} , 

where £>t denotes the elements of @ of order < t, as above. 

Proof. Let >- be the graded reverse lexicographic ordering in C[9] such that d\ -< 82 -< ■ • • -< d n . We 
consider the initial In(^) = {In(A) | A £ $)} of ffl for the monomial ordering It is a finite set of 
increasing monomials D := { d a °, d ai ,..., d 06 - 1 } ; which are the leading monomials of the elements 
of A = {Aq, Ai,..., A^_i} C S 1 . As 1 £ and is the lowest monomial >-, we have Aq = 1. As >- 
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is refining the total degree in C[0], we have ord(Aj) = |aj| and 0 = ord(A 0 ) < • • • < ord(A,s_i). 
Moreover, every element in @ t reduces to 0 by the elements in A. As only the elements Aj of 
order < t are involved in this reduction, we deduce that $>t is spanned by the elements Aj with 
ord(Aj) < t. 

Let E = {ckQ) ■ • •, & 5 - 1 }. The elements Aj are of the form 


A,- = 


aj! 




E ' 

I/3HKI 


«?■ 


By auto-reduction of the elements Aj, we can even suppose that /3 ^ E in the summation above, 
so that they are of the form d4|. 

Let S(£) ={(x-^)“°,...,(x-^) a ' 5 - 1 } c C[x]. As (Aj((x-£) a i))o<i,j<«-i is the identity matrix, 
we deduce that B is a basis of C[x]/Q, which is dual to A. 

As is stable by derivation, the leading term of ^(A j) is in D. If is not zero, then it 

is the leading term of ^-(A j), since the monomial ordering is compatible with the multiplication 
by a variable. This shows that D is stable by division by the variable di and that B is connected 
to 1. This ends the proof of the lemma. □ 


Such a basis of S> can be obtained from any other basis of S> by transforming first the coefficient 
matrix of the given dual basis into row echelon form and then reducing the elements above the 
pivot coefficients. The integration method described in [20] computes a primal-dual pair, such that 
the coefficient matrix has a block row-echelon form, each block being associated to an order. The 
computation of a basis as in Lemma 14.II can be then performed order by order. 

Example 4.2. Let 

fl = X! 

which has a multiplicity 3 root at £ = 
a primal-dual pair 

B = {1, xi,x 2 } , A = jl,<9i + d- 2 ,d 2 + E? + did 2 + ^ d\ 

This primal dual pair does not form an orthogonal pair, since ( d\ + d 2 )(x 2 ) 0. However, using 

let say the degree lexicographic ordering such that x\ > x 2 , we easily deduce the primal-dual pair 
of Lemma 14.11 

B = { l,xi,xf}, A = A. 

Throughout this section we assume that we are given a fixed primal basis B for C[x]/Q such 
that a dual basis A of satisfying the properties of Lemma f 1.1 1 exists. Note that such a primal 
basis B can be computed numerically from an approximation of £ and using a modification of the 
integration method of [201 . 

Given the primal basis B , a dual basis can be computed by Macaulay’s dialytic method which can 
be used to deflate the root £ as in [14] . This method would introduce n + (S — 1) (( n ^°) — $) new 
variables, which is not polynomial in o. Below, we give a construction of a polynomial system that 
only depends on at most n + nS(S — l)/2 variables. These variables correspond to the entries of 
the multiplication matrices that we define next. Let 

Mi : C[x]/Q -> C[x]/Q 

V ^ (Xi-£i)p 

be the multiplication operator by Xi — £j in C[x]/Q. Its transpose operator is 
M* : ->• @ 

A h> Ao Mi = (xi - £j) • A = E— (A) = d di s (A) 


- x 2 + x\,f 2 =xi- x 2 + x\, 

- (0,0). The integration method described in [20] computes 
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where Ld = Q^ C C[0]. The matrix of Mi in the basis B of C[x]/Q is denoted Mj. 

As B is a basis of C[x]/Q, we can identify the elements of C[x]/Q with the elements of the vector 
space span c (B). We define the normal form N(p) of a polynomial p in C[x] as the unique element 
b of span c (i?) such that p — bG Q. Hereafter, we are going to identify the elements of C[x]/Q with 
their normal form in span c (i3). 

For any polynomial p(x 1 ,..., x n ) G C[x], let p(M) be the operator of C[x]/Q obtained by replacing 
Xi - & by Mi. 

Lemma 4.3. For any p € C[x], the normal form of p is N(p) = p(M)(l) and we have 

P(M)(1) = A 0 (p) 1 + Ai(p) (x - 0 ai + • • ■ + A S -i(p) (x - O 04 - 1 . 

This shows that the coefficient vector [p] of N{p) in the basis B of is [p] = {Ai{p))o<i<g-i- 
The following lemma is also well known, but we include it here with proof: 

Lemma 4.4. Let B as in m and denote the exponents in B by E := {on,..., ag} as above. Let 

n 

E+ := IJt-E + e,) 

*=l 

with E + et = {( 71 ,..., ji + 1,... , 7 n ) : 7 G E} and we denote d(E) = E + \ E. The values of 
the coefficients v a ,p for (cc, ($) G E x d(E) appearing in the dual basis 0 uniquely determine the 
system of pairwise commuting multiplication matrices Mj, namely, for i = 1 ,... ,71 


0 

\ ,Gi 

^ OL2 

l/ a d _ 1 ,e i 

0 

0 

Vot2 

k'dd-i.ai+ei 

0 

0 

0 

v a d _i,OL d _ 2 +ei 

0 

0 

0 

0 


Moreover, 

_ I 1 if ai=a k + e, 

v 'a i ,a k +e j — S . 

[0 if a k + e 3 G E, ai a k + ej. 

Proof. As Ml acts as a derivation on and @ is closed under derivation, so the third property 
in Lemma 14.11 implies that the matrix of Ml in the basis of A of @ has an upper triangular form 
with zero (blocks) on the diagonal. 

For an element A ; of order k, its image by Ml is 

Mt{Aj) = ( Xi - fi) ■ Aj 

= £ Aj((xi — fi)(x — £) ai )Ai 

\ai\<k 

= J2 A,'((x- 0 “‘ +ei )A i = £ o aj , ai+ei Ai. 

\ai\<k 

This shows that the entries of M, are the coefficients of the dual basis elements corresponding to 
exponents in E x d(E). The second claim is clear from the definition of M,. □ 

The previous lemma shows that the dual basis uniquely defines the system of multiplication ma- 
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trices for i = 1 ,..., n 



> 

o 

-b) 

> 

Oo 

1 

IT 

I 

£rv 

> 

o 

((X- 


A 5 _! ((x-0“l+ e 

> 

o 

((x- 


A,i ((x-£) Q5+e 


0 

,e 

VOL2 


0 

0 

Va. 2 ,a:i +e i 

t/ a s _ 1 ,ai+e i 

0 

0 

0 

l / a s _ 1 ,a s _2+e i 

0 

0 

0 

0 


Note that these matrices are nilpotent by their upper triangular structure, and all 0 eigenvalues. 
As o is the maximal order of the elements of we have M 7 = 0 if |q| > o. 


Conversely, the system of multiplication matrices Mi,..., M„ uniquely defines the dual basis as 
follows. Consider v aijl for some (a,;, 7 ) such that |q| < o but 7 ^ E + . We can uniquely determine 
v aitl from the values of {v aj ,p '■ G E x d(E)} from the following identities: 

v ai)1 = A*((x - £) 7 ) = [M (x _| ) 7 ]i,i = [M 7 ] M . ( 6 ) 


The next definition defines the parametric multiplication matrices that we use in our constriction. 

Definition 4.5 (Parametric multiplication matrices). Let B as in Q, and E , d(E) as in Lemma 
14.41 We define array p of length nS(S— l )/2 consisting of 0’s, l’s and the variables p ai ,p as follows: 
for all a*, Ofc £ E and j £ {1,..., n} the corresponding entry is 


{ 1 if Oj = Ofe + ej 

0 if CXk + Gj £ E , on ^ Oik + Bj 

l^oti ,Qfc+ej if rrfc T Gj ^ E. 

The parametric multiplication 


nl(p) := 


matrices are denned for i = 1 ,..., n by 


0 

f^oti ,ei 

f^ot2,ei 

pas- 1 ,ei 

0 

0 

t^Ot2iOtl-\-Gi 

l^otd-i ,o;i+ei 

0 

0 

0 

fJ'Q'd-l,OL8-2+ei 

0 

0 

0 

0 


(7) 


( 8 ) 


We denote by 


M(/x) 7 :=M 1 (/i) 7 I ---M„(/x) 7 ”, 


and note that for general parameters values p, the matrices M, (/x) do not commute, so we fix their 
order by their indices in the above definition of M(/x) 7 . 


Remark 4.6. Note that we can reduce the number of free parameters in the parametric multipli¬ 
cation matrices by exploiting the commutation rules of the multiplication matrices corresponding 
to a given primal basis B. For example, consider the breadth one case, where we can assume that 
E = {0, ei, 2ei,..., (5 — l)ei}. In this case the only free parameters appear in the first columns of 
M 2 (/x), ..., M„(/x), the other columns are shifts of these. Thus, it is enough to introduce (n— 1)(A— 1) 
free parameters, similarly as in na. In Section [5] we present a modification of m Example 3.1] 
which has breadth two, but also uses at most (n — 1)(5 — 1 ) free parameters. 
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Definition 4.7 (Parametric normal form). Let K C C be a field. We define 
W z , m :K[x] -»■ K[z,/i]* 

P Af z ,„(p) : = E A 9 zb) M (M) 7 [l]- 

7SN n 

where [1] = [1,0,..., 0] is the coefficient vector of 1 in the basis B. This sum is finite since for 
M > S, M(/z ) 7 = 0, so the entries of W Zl/i (p) are polynomials in fi and z. 

Notice that this notation is not ambiguous, assuming that the matrices (/z) (i = 1 ,,n) are 
commuting. The specialization at (x, /z) = (£, u) is the vector 

•A k, v {p) = [A 0 (p), • • •, A 5 _i(p)]‘ g C 5 . 

4.1 The multiplicity structure equations of a singular point 

We can now characterize the multiplicity structure by polynomial equations. 

Theorem 4.8. Let K C C be any field, f G K[x] w and let f G C n be an isolated solution of f. Let 
M i(n) for i = 1 ,... n be the parametric multiplication matrices as in (01 and be the parametric 
normal form as in Defn. at z = £. Then the ideal J j of C[/z] generated by the polynomial 
system 

(A kAh) for k m l,... ,N, 

\m*( jLt) -Kjifi) -Mi(yu) ■ Mi (At) for — 

is the maximal ideal 

- v a ,p, (a,fi) £ Ex d(E)) 

where are the coefficients of the dual basis defined in 0 - 

Proof. As before, the system (0) has a solution n a< p = n a ,g for (a, fi) G E x d(E). Thus J £ C m„. 
Conversely, let C = C[/z]/ and consider the map 

$ : C[x]-► C*, p^fif^(p). 

Let K be its kernel. Since the matrices M,;(/z) are commuting modulo J^, we can see that K is an 
ideal. As fk £ K, we have T := (/*,) C K. 

Next we show that Q C K. By construction, for any a £ N n we have modulo 

A/U(x-£D= E V?((x-O q )M(m) 7 [ 1 ]=M(m)“[ 1 ]. 

Using the previous relation, we check that Vp,q £ C[x], 

®(pq) =p(£ + M(/i))$(g) (io) 

where p(£ + M(/z)) is obtained by replacing Xj — £, by M z (/u). Let q £ Q. As Q is the m^-primary 
component of I, there exists p £ C[x] such that p(f) / 0 and pq £ 1. By (ITT)]) . we have 

4>(p q) = p{f + M (/z))$(g) = 0 . 

Since p(£) ^ 0 and p(£ + M(/z)) = p{f)Id + N with N lower triangular and nilpotent, p(£ + M(/z)) 
is invertible. We deduce that $(( 7 ) = p(£ + M(/z)) _ 1 4>(pg) = 0 and q £ K. 

Let us show now that 4> is surjective and more precisely, that <E>((x — fi) ak ) = (abusing the 
notation as here has length S not n). Since B is connected to 1, either ctk = 0 or there exists 
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a.j £ E such that a k = ctj + e* for some i £ {1,... ,n}. Thus the j th column of Mj(/x) is e*, by 
0 . As {Mj(/x) : i = 1,... , n} are pairwise commuting, we have M(/x) afe = Mj(/x)M(/x)“ J ' , and if we 
assume by induction on \ctj\ that the first column of is ej, we obtain M(/x)“ fc [l] = e*,. Thus, 

for k = l,...,S, $((x-£)“ fc ) =e k . 

We can now prove that m„ C Jj. As Mi(y) is the multiplication by (xi — &) in C[x]/Q, for any 
b £ B and i = 1 ,n, we have (xj — fi)b = Mi(v)(b) + q with q £ Q C K. We deduce that for 
k = 0 ,..., 6 — 1 , 

$((*« - £*) (x - S)“ fc ) = Mi(/x)$((x - £)“*) = Mt*)(ek) = 

This shows that p a p — v a p £ for (a,/3) £ E x d(E) and that m„ = J^. □ 

In the proof of the next theorem we need to consider cases when the multiplication matrices do 
not commute. We introduce the following definition: 

Definition 4.9. Let IcCbe any field. Let C be the ideal of K[z,/x] generated by entries of the 
commutation relations: fh(/x) ■ M j(p) — M j(p) ■ Mj(/x) = 0, i, j = 1,... ,n. We call C the commutator 
ideal. 


Lemma 4.10. For any field K C C, p £ K[x], and i = 1 ,... ,n, we have 
A/' z ,/x(a up) = XiN^ip) + M i (/x)A/’ Z|#i (p) + O i:ll (p). 
where Oi^ : K[x] —> K[z, p] s is linear with image in the commutator ideal C. 
Proof. JV z ^{xip) = E 7 df(xip) M(/x) 7 [l] 

= E 4 3?(P) »W[ 1 ] + E ^ ST* ip) M(^) 7 [l] 

= X, E ^d2(p) M(^) 7 [l] + £ ±d2(p) M(/x) 7+Si [1] 


( 11 ) 


1 - 


= XiAf^ip) + M i(p) I E — d Zip) m (m) 7 [1] 


E ^ 5 z(p) °*.7(m)[1] 


where Ch j7 = Mi(/x)M(/x ) 7 — M(/x) 7+6i is a <5 x 6 matrix with coefficients in C. Therefore, : p H>• 

Z ) 7 ip) C>i ; 7 (/x) [1] is a linear functional of p with coefficients in C. 


□ 


The next theorem proves that the system defined as in m for general z has (£, v) as a simple root. 

Theorem 4.11. Let f £ K[x] jV and f £ C" be as above. Let M,;(/x) for i = 1,... n be the parametric 
multiplication matrices defined in (0 and JV X .^, be the parametric normal form as in Defn. \f?\ 
Then (z,/x) = (£, v) is an isolated root with multiplicity one of the polynomial system in K[z, p\: 


( 12 ) 


Wz, M (/fc) = 0 for k = 1,..., N, 

Mi(/x) ' m j (m) - Mj(m) ' Mi(/x) =0 for i, j = 1, .. . ,n. 

Proof. For simplicity, let us denote the (non-zero) polynomials appearing in (1121) by 

Pi,..., Pm GK[z,/x], 

where M < NS + n[n — 1)(<5 — 1)(<5 — 2)/4. To prove the theorem, it is sufficient to prove that 
the columns of the Jacobian matrix of the system [Pi,..., Pm] at (z, /x) = (£, v) are linearly 
independent. The columns of this Jacobian matrix correspond to the elements in C[z,/x]* 

di ,(,..., d n £, and d Ma for (a,/3) £ E x d(E), 
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where di^ defined in m for z replacing x, and d fla p is defined by 

( q ) = - 7 —^— | (z, M )=(l,s> for q £ C[z, p]. 

Suppose there exist ai,..., a n , and a ajj g S C for (a, (3) £ E x d{E ) not all zero such that 
A := ai9i^ + • • • + a n d n ^ + ( ia,pd ila p £ C[z, /j]* 

a, /3 

vanishes on all polynomials Pi,, Pm in (1121) . In particular, for an element Pi(p) corresponding 
to the commutation relations and any polynomial Q £ C[x, /i], using the product rule for the linear 
differential operator A we get 


A (PiQ) = A (Pi)Q(£, v) + Pi(v)A(Q) = 0 


since A (Pi) = 0 and Pi{v) = 0. By the linearity of A, for any polynomial C in the commutator 
ideal C, we have A (C) = 0. 

Furthermore, since A(A/’ Zi/i (/fc)) = 0 and 

A/f,^(/fc) = [A 0 (/fc), • ■ • .Aj_i {fk)Y, 


we get that 

(ai9i,e 4-h a n d n ,z) • A S -i(fk) + ^ p 7 (;z) dj(fk) = 0 (13) 

where p 1 £ C[/i] are some polynomials in the variables p that do not depend on fk- If a \,..., a n 
are not all zero, we have an element A of C[9j] of order strictly greater than ord(A, 5 _i) = o that 
vanishes on /i,..., /jv ■ 

Let us prove that this higher order differential also vanishes on all multiples of /*, for k = 1,..., N. 
Let p £ C[x] such that = 0, A(Af z ^{p)) = 0. By (fill) , we have 


■A4 tV {{xi - ii)p) 

= (*i - £,i)N^ u (p) + M i(u)Af iiV (p) + Oi, v (p) = 0 

and A(A4, M ((zi - &)p)) 

= A((xi - ii)M z ^{p)) + A(Mi(M)-A/z iAl (p)) + A(0 /i (p)) 

= A (xi - £,i)Afti,v{p) + (& - Ci)A(A/z iA1 (p)) 

+ A(M 4 (/x))A/^, m (p) + M, ; (r/)A(A/' z , A1 (p)) 

+ A (0,;, M (p)) 

= 0 . 


As N^, u (fk) = 0, A(A/’ z ,^i(/fc)) = 0, i = 1 ,..., N, we deduce by induction on the degree of the 
multipliers and by linearity that for any element / in the ideal / generated by fi,..., /at, we have 


A4^(/) = 0 and A(A/" z , m (/)) = 0, 

which yields As/- 1 . Thus we have As/- 1 - flC[c?^] = Q- 1 (by Lemma [All) . As there is no element 
of degree strictly bigger than o in Q ± , this implies that 


CL\ — * * * — d n — 0. 


Then, by specialization at x = £, A yields an element of the kernel of the Jacobian matrix of the 
system ©. By Theorem 14.81 this Jacobian has a zero-kernel, since it defines the simple point v. 
We deduce that A = 0 and (£, v) is an isolated and simple root of the system (fl^l) . □ 
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The following corollary applies the polynomial system defined in (11211 to refine the precision of an 
approximate multiple root together with the coefficients of its Macaulay dual basis. The advantage 
of using this, as opposed to using the Macaulay multiplicity matrix, is that the number of variables 
is much smaller, as was noted above. 

Corollary 4.12. Let f £ K[x] w and £ £ C" be as above, and let Ao(i'),..., be its dual 

basis as in W- Let E C N" be as above. Assume that we are given approximates for the singular 
roots and its inverse system as in 0 

£ = £ and u ai ,p — v ai4 3 Va* £ E, /3 & E, |/3| < o. 


Consider the overdetermined system in K[z, p] from 112 1) . Then a random square subsystem of \12\) 
will have a simple root at z = £, p = v with high probability. Thus, we can apply Newton's method 
for this square subsystem to refine £ and v ai ,p for (ccj ,/?) £ E x d(E). For v aitl with 7 ^ E + we 
can use (01 for the update. 


Example 4.13. Reconsider the setup from Ex. l3.21 with primal basis {1, X 2 } and E = {(0, 0), (0,1)}. 
We obtain 


Mt) = 


0 0 
p 0 


and M 2 (/x) = 


0 0 
1 0 


The resulting deflated system in (fl2l) is 


F(Z!,Z2,p) 


Z\ + Z 2 
P + 2Z2 

2 1 2 

Zi + Z2 

2pz\ + 2z 2 


which has a nonsingular root at ( 27 , Z 2 , p) = ( 0 , 0 , 0 ) corresponding to the origin with multiplicity 
structure { 1 , 82 }. 


5 Examples 

Computations for the following examples, as well as several other systems, along with Matlab 
code can be found at www.nd.edu/~jhauenst/deflation/ 

5.1 A family of examples 

We first consider a modification of E23 Example 3.1]. For any n > 2, the following system has n 
polynomials, each of degree at most 3, in n variables: 

3,223,2 3,2 2 

*^1 T 37 X2 “b X2 £3, - - - , %n —1 ~b *^n—1 % n . 

The origin is a multiplicity 5 := 2 n root having breadth 2 (i.e., the corank of Jacobian at the origin 
is 2). 

We apply our parametric normal form method described in § 0 Similarly as in Remark 14.61 we 
can reduce the number of free parameters to be at most [n — 1)(<5 — 1 ) using the structure of the 
primal basis B = {x^x^ : a < 2 n_1 , b < 2}. 

The following table shows the multiplicity, number of variables and polynomials in the deflated 
system, and the time (in seconds) it took to compute this system (on a iMac, 3.4 GHz Intel Core i7 
processor, 8 GB 1600Mhz DDR3 memory). Note that when comparing our method to an approach 
using the null spaces of Macaulay multiplicity matrices (see for example [B, .I3])j we found that for 
n > 4 the deflated system derived from the Macaulay multiplicity matrix was too large to compute. 
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so the size of the Macaulay multiplicity matrix 


This is because the nil-index at the origin is 


is n • 


+n—1 
n —1 


) x r;n 


2 n "\ 



New approach 

Null space 

n 

mult 

vars 

poly 

time 

vars 

poly 

time 

2 

4 

5 

9 

1.476 

8 

17 

2.157 

3 

8 

17 

31 

5.596 

192 

241 

208 

4 

16 

49 

100 

19.698 

7189 

19804 

> 76000 

5 

32 

129 

296 

73.168 

N/A 

N/A 

N/A 

6 

64 

321 

819 

659.59 

N/A 

N/A 

N/A 


5.2 Caprasse system 

We consider the Caprasse system mm- 

4xi 2 + 

4xiX3 + 

at the multiplicity 4 root £ = (2, —y/—3, 2, s/—3). 

We first consider simply deflating the root. Using the approaches of la ni na , one iteration 
suffices. For example, using an extrinsic and intrinsic version of EUH, the resulting system 
consists of 10 and 8 polynomials, respectively, and 8 and 6 variables, respectively. Following flO] , 
using all minors results in a system of 20 polynomials in 4 variables which can be reduced to a 
system of 8 polynomials in 4 variables using the 3x3 minors containing a full rank 2x2 submatrix. 
The approach of §[3] using an |i| = 1 step creates a deflated system consisting of 6 polynomials in 
4 variables. In fact, since the null space of the Jacobian at the root is 2 dimensional, adding two 
polynomials is necessary and sufficient. 

Next, we consider the computation of both the point and multiplicity structure. Using an intrinsic 
null space approach via a second order Macaulay matrix, the resulting system consists of 64 
polynomials in 37 variables. In comparison, using the primal basis {l,Xi,X 2 , X 1 X 2 }, the approach 
of § |4] constructs a system of 30 polynomials in 19 variables. 

5.3 Examples with multiple iterations 

In our last set of examples, we consider simply deflating a root of the last three systems from El 
§ 7] and a system from T2J § 1], each of which required more than one iteration to deflate. These 
four systems and corresponding points are: 

1: (xf — X 2 X 3 X 4 , x| — X 1 X 3 X 4 , x| — X 1 X 2 X 4 , x\ — X 1 X 2 X 3 } at (0, 0, 0,0) with 5 = 131 and o = 10; 

2: {x 4 , x 2 y + y 4 ,z + z 2 — 7x 3 — 8 x 2 } at (0, 0, —1) with 5 = 16 and o — 7; 

3: {14x + 33y — 3\/5(x 2 + 4 xy + 4 y 2 + 2) + s/7 + x 3 + 6 x 2 y + 12 xy 2 + 8 y 3 , 41x — 18 y — -s/5 + 8x 3 — 
12x 2 y + 6 xy 2 — y 3 + 3s/7(4xj/ — 4x 2 — y 2 — 2)} at Z 3 « (1.5055, 0.36528) with <5 = 5 and o = 4; 

4: {2xi + 2x 2 -f- 2x2 + 2 x 2 + x 2 — 1, (xi + X 2 — X 3 — l ) 3 — x 3 , 

(2x? + 5 x 2 + IOX 3 + 5x§ + 5 ) 3 — lOOOxf} at (0, 0, —1) with <5 = 18 and o = 7. 

We compare using the following four methods: (A) intrinsic slicing version of [6lll.3j: (B) isosingular 
deflation [TO] via a maximal rank submatrix; (C) “kerneling” method in [8j; (D) approach of 
§ [3] using an |i| = 1 step. We performed these methods without the use of preprocessing and 


/(x i,x 2 ,x 3 ,x 4 ) = 


Xl X 3 — 4 xiX 2 2 X 3 — 4 xi 2 X’ 2 X 4 — 2 X 2 3 X 4 

10 X2 2 — 4 X1X3 + 10 X2X4 — 2 , 

X1X3 3 — 4X2X3 2 X4 — 4X1X3X4 2 — 2X2X4 3 

IOX2X4 — 4 x 3 2 + 10 X4 2 — 2 , 
X 2 2 %3 + 2X1X2X4 — 2 xi — X3, 

X4 2 Xi + 2 X2X3X4 — 2 X3 — Xi 
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postprocessing as mentioned in § [3] to directly compare the number of nonzero distinct polynomials, 
variables, and iterations for each of these four deflation methods. 



Method A 

Method B 

Method C 

Method D 


Poly 

Var 

It 

Poly 

Var 

It 

Poly 

Var 

It 

Poly 

Var 

It 

1 

16 

4 

2 

22 

4 

2 

22 

4 

2 

16 

4 

2 

2 

24 

11 

3 

11 

3 

2 

12 

3 

2 

12 

3 

3 

3 

32 

17 

4 

6 

2 

4 

6 

2 

4 

6 

2 

4 

4 

96 

41 

5 

54 

3 

5 

54 

3 

5 

22 

3 

5 


For breadth one singular points as in system 3, methods B, C, and D yield the same deflated system. 
Except for methods B and C on the second system, all four methods required the same number of 
iterations to deflate the root. For the first and third systems, our new approach matched the best 
of the other methods and resulted in a significantly smaller deflated system for the last one. 
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